Quantum scattering model of energy transfer in photosynthetic complexes 
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We develop a quantum scattering model to describe the exciton transport through the Fenna- 
Matthews-Olson(FMO) complex. It is found that the exciton transport involved the optimal quan- 
tum coherence is more efficient than that involved classical behavior alone. Furthermore, we also 
find that the quantum resonance condition is easier to be fulfilled in multiple pathways than that 
in one pathway. We then definitely demonstrate that the optimal distribution of the pigments, the 
multitude of energy delivery pathways and the quantum effects, are combined together to contribute 
to the perfect energy transport in the FMO complex. 

PACS numbers: 87.15.A-; 71.35.-y; 87. 15. hj 



Photosynthesis is a fundamental biological process that 
harvests solar energy to power life on Earth. It starts 
with the absorption of a photon of sunlight by the an- 
tenna molecules, followed by transfer of the energy to the 
reaction centers. Remarkably, the energy transfer from 
the antennas to the reaction center is almost perfectly 
efficient fl. Then an intriguing question arises: what 
makes the energy transfer so efficient? Up to date, there 
are two acceptable contributing factors [2j: The first one 
is optimal space distribution of the pigments in light- 
harvesting complexes [3|, where the optimal space is just 
closely enough to enable fast energy transfer, but far 
enough apart to prevent the molecular orbital of the 
pigments from overlapping. The second factor is the 
supramolecular organization of photosynthetic appara- 
tus, which allows the multiple energy delivery pathways 
to connect to the reaction centre 

Recently, the experimental and theoretical studies 
show that there exists long-lived quantum coherence in 
this energy transfer in several photosynthetic complexes 
, which suggests that quantum coherence could be a 
third factor in optimizing energy transfer efficiency. How- 
ever, to what extent quantum effects contribute to the 
efficiency of energy transfer is largely unknown. There- 
fore, it would be significant important to establish quan- 
titatively whether or not photosynthetic light harvesting 
that involves quantum coherence is truly more efficient 
than it would be using classical mechanisms alone Q. 

In this paper, we develop a quantum scattering model 
to describe the exciton (the energy carrier) transport in 
the FMO complex. We explore respectively the quan- 
tum and classical transports of the exciton in one, two 
and multiple pathways and show that the quantum trans- 
port under the resonant condition can lead to perfect 
efficiency. Since the resonance is due to the quantum 
coherence and it is impossible in classical physics, we 
thus demonstrate definitely that the exciton transport in- 
volved optimal quantum coherence is more efficient than 
that for the classical case. We show that the resonant 
transport with perfect efficiency occurs when the quan- 



tum phase of the exciton accumulated in the pathways 
are fulfilled with the resonant condition. Since the phase 
of the exciton is determined by the space distribution of 
the pigments, the optimal distribution of the pigments 
is then indeed significant in the near-unity transfer ef- 
ficiency. Although the perfect efficiency is also possible 
for one pathway alone, the multiple pathways are essen- 
tial because the resonant region in multiple pathways is 
much larger than that in one pathway. The large reso- 
nant region is the mechanism that the excitons with a 
broad range of the spectrum ( electron absorbed differ- 
ent colors of light) can transport with perfect efficiency. 
Therefore, the three mentioned contributing factors, the 
optimal distribution of the pigments, the multitude of 
energy delivery pathways and the quantum effects, are 
combined together in our proposed model to contribute 
the perfect energy transport in FMO. 




FIG. 1. (Color online) The schematic representation of the 
light- harvesting apparatus of green sulphur bacteria, (a) The 
light is absorbed by the chlorosome antenna and then trans- 
ported to the reaction center through the FMO complex, 
(b) Each monomer has seven pigments labelled by ® — (7). 
Two main energy transfer pathways La and Lb are indi- 
cated by the red and green arrows 11], La- ®— 
Lb' a ring structure contains two paths @— ^(D— ^@ and 

Quantum scattering model of the FMO complex. The 
architecture of antenna light-harvesting complexes varies 
widely among photosynthetic organisms. A well-studied 
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example is the water-soluble FMO complex of green sul- 
fur bacteria. As shown in Fig. 1, the FMO complex 
comprises three identical monomers that each contains 
seven pigments. The seven pigments (Fig. 1(b)) can 
be combined into two major energy transport pathways 
[Tl| . A light pulse absorbed by the chlorosome antenna 
could excite an electron into a state called an exciton 
which is a quasi-particle and can be considered as the 
energy carrier in light- harvesting complexes. The exci- 
ton then transports to the reaction centre through FMO 
with nearly perfect efficiency. 

Exciton transport is accomplished in FMO through 
tunneling, where the exciton jumps from one pigment to 
another. In the previous theories (ji [ToL fl2| . each pigment 
is usually modeled by a two-level molecule and then the 
FMO complex are described by some coupled two-level 
molecules interacting with their environmental phonons. 
In this paper weproposed a different model, the scat- 
tering model [HI, Q| (see a brief review in Supplemen- 
tary material), to explain transmission efficiency of FMO. 
There are two basic elements in the scattering model, the 
scatterer and the ideal channel; they are connected with 
certain orders to form a quantum network. We denote 
the transfer matrix for the nth scatterer as M n and the 
transfer matrix for the ideal channel with the length I as 
XJ\. They can be written as 



Ui = 



ikl 





a ikl 



(i) 



where t n =T n 2 e llfn , r n =—iRne llfn 1 T n is the transmission 
probability of the nth scatterer, R n =l-T n is its reflection 
probability, <p n is the phase change in the transmitted 
wave and k is the wave number of the exciton. 

The transport of FMO can be well explored with the 
help of the transfer matrices in Eq. (1) as well as the 
junctions S and G to be defined below. Each pigment 
can be thought of as a potential well and there exists a 
potential barrier between two nearest neighbor pigments 
along the pathways. The potential barrier can be well 
described by the transfer matrix of the scatterer defined 
in Eq.(pQ). The wave function of the exciton in pigment 
j could be written as ipj(x)=Cje' lkx -\-Dje~ lkx , where Cj 
and Dj are the amplitudes of the wave function. The 
amplitudes Cj and Dj are related to its nearest neighbor 
pigment through the transfer matrix of the scatterer be- 
tween the two pigments. In addition, the transfer matrix 
of the pigment is the one as the ideal channel defined in 
Eq.(pQ). Therefore, we can describe the exciton transmis- 
sion through the whole FMO complex with the following 
three steps: (i) the two major energy transport pathways 
in each monomer is described by two different scattering 
pathways La and Lb, as shown in Fig. 2(a). (ii) Each 
monomer is represented by a ring with two branches La 
and Lb and connected to two leads through two junctions 
S. The transfer matrix of the monomer j (j=l,2,3) is 
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FIG. 2. The quantum scattering model we proposed to de- 
scribe the excitation transmission of FMO. (a) The scattering 
matrix representation of two major pathways La and Lb- (b) 
Each monomer is represented by the symbol Fj made by a ring 
with two branches La and Lb and connected to two leads, (c) 
The scattering matrix representation of the whole FMO com- 
plex. The symbols near the arrows denote the corresponding 
transmitted and reflected amplitudes. 



denoted as Fj, as show in Fig. 2(b). (hi) The transport 
of the whole FMO can be described in Fig. 2(c), where 
three monomers represented by Fj are connected to the 
antennas and reaction center through two junctions G. 
The exciton created at the antennas injects into FMO 
with 1 = 1 through the G junction connected with the 
antennas. For simplicity, we assume that there is no re- 
flection from the reaction center (the amplitude O in Fig. 
2(c) vanishes); In this case, the transmission coefficient 
of the whole FMO complex is then determined by the 
probability \0 ,2 \ indicated in Fig. 2(c). 

Transfer matrices of the FMO complex. We now turn 
to study the transfer matrices of the basic elements La, 
Lb, Fj and the whole FMO complex. The pathway La is 
represented by two scatterers connected one ideal channel 
and its transfer matrix is then given by 



L A = M 2 U i M 1 



(2) 



where I is the length of the 2nd pigment. The pathway 
Lb is composed of three branches connected by two S 
junctions (© and @). For each branch, we have 



pR 

rn 



pL 

m 



,(m = 1,2,3) 



(3) 



where T^Ui.M^MsUi^ T 2 =Ui 2 M 6 Ui 2 M 5 Ui 2 and 
Ts=MyUi 3 . Here l m is the distance between the scatter- 
ers or the distance between the scatterers and the junc- 
tions for the rath branch. The amplitudes of the three 
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outgoing waves of the junction S are related by a S ma- 
trix to the amplitudes of the three incoming waves, i.e., 



p? 

pi 




p 3 L 
of 



S\P?\, (4) 



where S is a 3 x 3 unitary matrix and its detailed form 
is given by Eq. (14) in Supplementary material. 

From Eqs. (1-4), the amplitudes C£,£>£ (a=A,B) of 
the wave function on the left of the pathway and the 
amplitudes C^D^ on the right of the pathway can be 
related by the transfer matrixes L ai that is, 



r I c <* 

L « I £)L 



(5) 



Therefore, instead of an array of the scatterers, each 
pathway may be replaced by an effective transfer matrix. 

We now study the exciton transport in one monomer of 
the trimer. The schematic representation of the monomer 
is shown in Fig. 2(b). Two pathways La and Lb form 
a ring and connected to two leads through two S junc- 
tions. With the help of the effective transfer matrix for 
the pathway a of the monomer j, we can obtain 



C 



D 



D 



r 



(7 = 1,2,3) (6) 



where Lj^ a =Ui ja Mj^ a Ui ja , and is the distance be- 
tween the scatterers and S junctions for the pathway a 
of the monomer j. 

The amplitudes of the three outgoing waves are related 
to the amplitudes of the three incoming waves through 
S matrix given by 




= S 




y 

= s I a 



i 



(7) 

From Eqs. (6-7), the amplitudes of the wave function 
on the left of the monomer j and the amplitudes on the 
right are related by the matrix Fj given by 



Vi 



i 



Vi 



(8) 



Similar to the transmission of one monomer, we can 
describe the exciton transport through FMO by two G 
junctions as shown in Fig. 2 (c). The amplitudes of the 
four outgoing waves of the junction G are related by a G 
matrix to the amplitudes of the four incoming waves 



x i 
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where G is a 4 x 4 unitary matrix given by Eq. (20) in 
Supplementary material. By solving Eq.® with 1=1 



and O=0, one can obtain the transmission coefficient 
T FM o = \0' 2 \ of the FMO complex. 

To understand the features of the transmission, we 
compare the transmission coefficients for different num- 
ber of pathways. For simplicity, we assume that all scat- 
terers have the same transmission probability T n =T and 
the phase change in transmitted wave (p n =0. Further- 
more, we compare the transmission difference between 
the classical and quantum cases. As for the classical 
transmission, it means that the exciton is a classical par- 
ticle without wave-like phase and thus the particle prob- 
ability is determined by the addition of the probability 
from different pathways and not the amplitudes of the 
wave function; In addition, the scatterer is determined 
uniquely by the transmission probability T. 

Transmission of the pathways La and Lb- Based on 
Eqs. (1,2,5), we can obtain the total transmission prob- 
ability Ta=\Ga\ 2 for pathway La under the conditions 
C\=\ and D%=0, 



T A = T 2 /[T 2 +4(1-T) cos 2 0a] 



(10) 



where (f>A=kL It is notable that the resonant tunneling 
Ta=1 appears under the condition (j) A =(n + 1/2)tt with 
n being an integer. However, as for the classical case, 
after considering the multi-reflection between the scat- 
terers Mi and we may obtain the total transmission 
probability T^=j^. 

Figure 3(a) shows the total transmission probability as 
a function of the path phase for both quantum and clas- 
sical case at strong scattering limit T=0.01. For classical 
case, T%=5 x 10 -3 , it is impossible for the particle to pass 
through the system at strong scattering limit. However, 
for quantum case, there exist two resonant peaks at which 
the total transmission probability takes its maximal value 
Ta=1. Therefore, the transport involved optimal quan- 
tum coherence is more efficient than that involved the 
classical behaviors alone. 

For pathway Lb, we set Cjj = l, Z}^=0, kli=kl2=4>B 
and kls=(j) B . From Eqs. (2-5), the total transmission 
probability Tb=\C§\ 2 is given by 



T B = 



T 2 sin 2 7+ [(1 + R) cos7 + 2i? cos D r 



(11) 



where 7=20^ The resonant tunneling 
Tb=1 will happen under the conditions 7=7r and 
<\) B = \ arcsin[(ri+T-2)/2i?]. 

The path phases 0a, <\>b and <p B denote the dis- 
tance between the neighbor scatterers or the distance 
between the scatterers and junctions, which are deter- 
mined actually by the space distribution of the seven 
pigments. When the pigments are optimally spaced 
such that (j)A=\^ for La and 7=7r, <Pb = \ arcsin[(T2 +T- 
2)/2R] for Lb, the quantum resonant tunneling will 
happen and the exciton can pass through the barriers 
with perfect efficiency. Therefore, we have theoretically 
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FIG. 3. (Color online) The transmission probability as a 
function of the path phases at T=0.01. (a) Single path- 
way La for both quantum and classical cases. (b) A 
monomer with double pathways, (c) The whole FMO with 
multiple pathways at ^s^^arcsin (V0A2T-VR). (d) The 
resonant structures: (I) one pathway La, (II) two path- 
ways at (\)b—\ arcsin(v / L-V r R), (III) multiple pathways at 
0i+02=tt/2 and 3 = ^ arcsin(V0.12T-^R). 



demonstrated that optimally spaced pigments in light- 
harvesting complexes can facilitate the energy transfer. 

Transmission of one monomer. From Eqs. (6-8), 
we can obtain the total transmission probability of one 

and 



r,R\2 



monomer Tm- 
klj,B=4>B, which is given by 



at x j - 



1, 2/f =0, kl jyA ~- 



Tm 



H 2 /(P 2 + Q 2 ) 



(12) 



where H=VT(2VR+sm 20A+sin 2<j> B ),P=2 sin 2 
-sin 2 (^A-^B)+V^(sin20A+sin2(/)B),Q=sin2((/)A+0 J B) + 
V r R(cos20A+cos2^B)- We can find that the resonant 
tunneling Tm=1 under the conditions 0a +</>#=§ and 
4>a = \ arcsin(v / T-V^R). In contrast, the classical trans- 
mission probability is less than its upper bound 

3T 



1+2T' 



Figure 3(b) shows the total transmission probability 
of one monomer versus the path phases at T=0.01. For 
classical case, Tc<0.029, the particle passes through the 
monomer with very small probability. However, for quan- 
tum case, there are eight resonant regions where the par- 
ticle can pass through the monomer with perfect effi- 
ciency. Compared with single pathway, the resonant re- 
gions for two pathways become larger and then it is much 
easier for the exciton absorbed different colors of light to 
pass through the monomer. 

Transmission of the whole FMO complex. From Eqs. 
(6-9), we can obtain the total transmission probability 
Tfmo = \0 1 2 of the whole FMO complex. For simplicity, 
we assume klj y A+klj,B = ^ and the phases accumulated 



between G and S are an integer of 2tt. Then we have 

T FM o = H 2 /Q 2 , (13) 

where H=6VT qi q 2 q 3 (3R + 2p x ^R + \p\ - ±p 2 ), 
Q = -6 Pl R? -(Hp? - 6p 2 + 9)R 2 - 2[30p - 
7ps + pi (6 + 7p 2 )]R% - {lp\ + 34p Pi + 4p 2 - 
3p 2 - 4p 4 )# + (9p + 2pi)(p 2 - p\)R^ - 9p 2 - 
2piPo - \v\ + \pa- Here p =sin2(/>i sin2(/> 2 sin2(/> 3 , 
p n =sin n 2(/>i+sin n 20 2 +sin n 20 3 (n=l, 2, 3, 4), and 
qj=sin2(j)j+\fR with (f)j=hlj^A- It is found that 
the exciton can resonantly pass through FMO when 
01=02=03=^ arcsin(v / T-V r R). In contrast, the total 
classical transmission probability Tp MO is less than its 
upper bound ■ Figure 3(c) shows the quantum 

transmission of the whole FMO complex for fixing 
03=^ arcsin(V0.12T-V^) and T=0.01. Compared with 
the transmission of double pathways shown in Fig. 3(b), 
the resonant region of the whole FMO becomes larger. 
However, the transmission Tp MO <0.083 for classical 
case, which is pretty small. 

Before ending, we make several remarks on the main 
results that can be obtained from the scattering model: 
(i) The mechanism underlying the perfect energy effi- 
ciency is the resonant transmission induced by the quan- 
tum coherent effects, which is impossible for the classi- 
cal mechanism, (ii) Multi-pathway plays significant role 
in achieving remarkable efficiency of the energy trans- 
fer. To understand this point, we plot in Fig. 3(d) the 
transmission coefficients for one, two, and multiple path- 
ways. The resonant regions for FMO with multiple path- 
ways almost occupy all parameter space which ensures 
that the exciton can pass through FMO with perfect effi- 
ciency. The large resonant region means that the exciton 
absorbed different colors of light can still pass through 
FMO with high efficiency, (iii) The resonant conditions 
are determined by the space distribution of the pigments. 
When the pigments are optimally spaced such that the 
quantum resonant tunneling will happen, the exciton can 
pass through FMO with perfect efficiency. Therefore, 
our results prove the conjecture that quantum coherence 
might 'wire' together the final energy acceptors in an- 
tennas, thus compensating for the relatively weak cou- 
pling between the pigments in FMO [2|, lS|. In addition, 
FMO with multiple pathways has fault-tolerant function. 
When one pathway or one subunit is blocked, the exciton 
can still pass through other pathways. 

In summary, we have proposed a scattering model to 
compare the quantum and classical energy transfer in 
FMO. We found that, the aforementioned three factors, 
the optimal space distribution of the pigments, multiple 
pathways, and quantum effects are combined together to 
ensure high efficiency of the energy transfer in FMO. 
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